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ABSTRACT 

The central regions of many interacting and early-type spiral galaxies are actively forming 
stars. This process affects the physical and chemical properties of the local interstellar medium 
as well as the evolution of the galaxies. We observed near-infrared H2 emission lines: v — 
1-0 S{1), 3-2 S{3), 1-0 5(0), and 2-1 S'(l) from the central ^ 1 kpc regions of the archetypical 
starburst galaxies, M82 and NGC 253, and the less dramatic but still vigorously star-forming 
galaxies, NGC 6946 and IC 342. Like the far-infrared continuum luminosity, the near-infrared 
H2 emission luminosity can directly trace the amount of star formation activity because the 
H2 emission lines arise from the interaction between hot and young stars and nearby neutral 
clouds. The observed H2 line ratios show that both thermal and non-thermal excittion are 
responsible for the emission lines, but that the great majority of the near-infrared H2 line emission 
in these galaxies arises from energy states excited by ultraviolet fluorescence. The derived physical 
conditions, e.g., far- ultraviolet radiation field and gas density, from [Cii] and [Oi] lines and far- 
infrared continuum observations when used as inputs to photodissociation models, also explain the 
luminosity of the observed H2 1-0 S{1) line. The ratio of the H2 1-0 5(1) line to far-IR continuum 
luminosity is remarkably constant over a broad range of galaxy luminosities; Lh2/ Lfib. — 10^'^, 
in normal late- type galaxies (including the Galactic center), in nearby starburst galaxies, and in 
luminous IR galaxies (LIRGs: Lpm > 10"'^"'^ i©)- Examining this constant ratio in the context 
of photodissociation region models, we conclude that it implies that the strength of the incident 
UV field on typical molecular clouds follows the gas density at the cloud surface. 

Subject headings: galaxies: individual (M82, NGC 253, NGC 6946, IC 342) — galaxies: ISM — galaxies: 
spiral — galaxies: starburst — infrared: ISM — ISM: lines and bands 
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1. INTRODUCTION 

The Infrared Astronomical Satellite (IRAS) de- 
tected a significant population of galaxies that 
emit large fractions of their bolomctric energy in 
the mid- and far- infrared (IR) bands. The ma- 
jority of luminous IR galaxies (LIRGs: Lfir > 
10^^ Lq) are powered by massive stars that have 
recently formed in their central kiloparsec (Gen- 
zel et al. 1998). The far-IR luminosities of these 
galaxies therefore measure the amount of star for- 
mation activity (Kennicutt 1998). Local starburst 
galaxies with far-IR luminosities of 10^" — 10^^ 
Lq and normal spiral galaxies with Lpiu — 10^ 
Lq also show evidence for active star formation in 
their central (< 1 kpc) regions, even though their 
far-IR luminosities are much lower than those of 
LIRGs (Thronson et al. 1987; Stacey et al. 1991). 

The formation of O and B stars affects the phys- 
ical and chemical state of the interstellar medium 
through stellar UV radiation and through the 
effects of stellar winds and supernovae (Telesco 
1988). Strong near-IR H2 emission lines have been 
observed from the central regions of galaxies, es- 
pecially from starburst galaxies and LIRGs. Since 
H2 vibration-rotation emission can trace the in- 
teraction of the neutral ISM with the UV radia- 
tion from hot young stars and its interaction with 
sources of dynamical energy such as supernovae, 
the H2 luminosity may also be a good tracer of 
star formation activity. The main motivation of 
this paper is to explore how well we can use this 
additional, readily observable tracer of the inter- 
action of massive stars with the dense interstellar 
medium to study massive star formation in the 
inner parts of galactic disks. 

The observational part of this paper emphasizes 
new, large-beam observations of emission in the 
vibration-rotation transitions of H2 in the near- 
IR from nearby late-type galaxies and low lumi- 
nosity starbursts. For normal galaxies, observa- 
tions of near-IR H2 lines have been restricted to 
local bright regions because the average surface 
brightness of the H2 emission is low. In a previ- 
ous study, we used an instrument optimized for 
observations of low surface brightness line emis- 
sion in the near-IR to detect enormously extended 
3° or ~ 400 pc) H2 emission along the plane of 
the Galaxy near the Galactic center (Pak, Jaffe, 
& Keller 1996a,b). The integrated and dered- 



dened H2 v — l OS'(l) luminosity from the ex- 
tended region (8 x 10"^ Lq) dominates over that 
from the previously known 5 pc circumnuclear gas 
ring (4.5 x 10^ Lq, Gatley et al. 1984). 

Our results for the inner 400 pc of the Milky 
Way motivated us to begin a study of H2 emission 
from the central kiloparsec of other nearby galax- 
ies where it may be possible to detect extended, 
low surface brightness H2 emission using the same 
large- aperture Fabry-Perot spectrometer used in 
the Galactic work. The large-beam observations 
we present here make it possible to compare for 
the first time the H2 emission to the far-IR con- 
tinuum and far-IR fine-structure line emission at 
similar angular resolution and over similar angu- 
lar scales. The observations also give us measure- 
ments of H2 line strengths from nearby objects 
with linear scales comparable to those obtained 
with slit spectrometers for the more distant lu- 
minous IR galaxies. With this broad range of 
spatially compatible observations, we investigate 
whether it is possible to formulate a generalized 
star formation law which can be applied to galax- 
ies over a wide range of levels in star formation 
activity. 

In order to analyze the H2 line data and to 
compare them with other star formation tracers, 
we need to identify the excitation mechanism re- 
sponsible for the H2 emission. In the central re- 
gions of galaxies, however, the excitation mecha- 
nism for the H2 emission has long been the subject 
of intense controversy (see the reviews of Mouri 
1994; Goldader et al. 1997a). Puxley, Hawarden, 
& Mountain (1988) observed nearby spiral galax- 
ies using a circular variable filter (A/AA = 120) 
and a beam full- width at half maximmn (FWHM) 
of 19'.'6. They detected 1 5'(1), 2 1 3(1), and 
1 S{0) lines and concluded that the H2 line pop- 
ulations are non-thermal. However, Moorwood 
& Oliva (1990) presented high spectral resolution 
data (A/AA — 1500) and pointed out that the 
2-1 S{1) flux measurement of Puxley, Hawarden, 
& Mountain (1988) could be affected by unre- 
solved stellar absorption lines. They argued that 
the H2 emission is mostly from shocks. Mouri 
(1994) suggested that high spectral resolution ob- 
servations are essential to discriminate between 
the excitation mechanisms. 

The ultimate energy source responsible for the 
excitation of the H2 vibration-rotation emission 
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can lie cither in the absorption of UV photons 
from the ambient interstellar radiation field or 
from local OB stars, or in shock waves excited 
by supernovae, protostellar flows, or other ener- 
getic events (Black & van Dishoeck 1987; Burton 
1992). Hydrogen molecules that absorb far-UV 
radiation make a transition to an excited elec- 
tronic state. At low densities, this transition is 
followed by fluorescence to the ground electronic 
state and a cascade through that state's vibration- 
rotation levels (Sternberg 1989; Pak et al. 2003). 
This interaction between H2 molecules and UV 
radiation takes place at the edges of molecular 
clouds where hydrogen makes the transition from 
predominantly atomic to predominantly molecular 
gas, in the photodissociation region (PDR: Tielens 
& Hollenbach 1985; Hollenbach & Tielens 1999). 
If uh < 5 X 10^ cm~^, the relative intensities of 
vibration-rotation transitions in the near-IR aris- 
ing in UV-excited H2 are insensitive to density or 
to UV field strength (Black & van Dishoeck 1987) 
and depend only on the branching ratios in the 
downward cascade. In dense {uh > 5 x 10^ cm~"^) 
gas heated by intense far-UV radiation, Sternberg 
& Dalgarno (1989) showed that inelastic H— H2 
and H2— H2 collisions can heat the gas at the cloud 
edge to temperatures approaching 1000 K. A com- 
bination of the coUisional de-excitation and col- 
lisional excitation in the high temperature layer 
alters the H2 level populations in this UV-excited 
gas. The expected changes in relative line strength 
with density have been verified by observations of 
molecular cloud edges (Luhman et al. 1997). In 
shocked regions, the H2 vibration-rotation energy 
level populations in the ground electronic state are 
thermalized by collisions and the level populations 
are usually in LTE (Shull & Beckwith 1982; Black 
& van Dishoeck 1987). Because the level popula- 
tions are thermalized and reflect a similar range 
of temperatures, shock- heated gas and hot, dense 
gas heated by UV radiation are difficult to distin- 
guish. 

Most previous observations of the central re- 
gions of galaxies showed that the H2 line ratios 
of 2-1 5(1) to 1-0 5(1) lie between the ratios 
one would expect for fluorescently excited H2 and 
what one would see in observations of UV-heated 
but collisionally excited gas or of gas collision- 
ally heated in shocks (Doyon, Wright, & Joseph 
1994; Kawaraet al. 1990; Koornneef & Israel 1996; 



Lumsden, Puxley, & Doherty 1994; Moorwood 
& Oliva 1990; Puxley, Hawarden, & Mountain 
1988, 1990; Davies et al. 2003). These intermedi- 
ate ratios suggest that both pure UV-fluorescence 
and thermalized UV-excited (or shock excited) gas 
are responsible for the H2 emission (Davies et al. 
2003). Detection of H2 lines arising from upper vi- 
brational states with ti > 2 is a strong indication 
of the presence of UV-fluorescent excitation since 
the high energies of the upper states make it very 
hard to excite these transitions thermally. There 
are marginal detections of such transitions from a 
few galaxies, e.g., a giant Hii region NGC 5451 
in MlOl (Puxley, Ramsay Howat, & Mountain 
2000), Arp 299 (Sugai et al. 1999), and some ul- 
trahiminous IR galaxies (Davies et al. 2003). An- 
other method for discriminating between the H2 
excitation mechanisms is to measure the H2 nu- 
clear spin degeneracy ratio (ortho- to para-ratio), 
which may deviate from its formation value if the 
gas is thermally processed, but remains constant 
in the case of non-thermal excitation. Observa- 
tions of 1-0 5(1) and 1-0 5(0) showed that the H2 
emission in the galactic centers arises mostly from 
UV-fluorescence (Harrison et al. 1998; Sugai et al. 
1999; Puxley, Ramsay Howat, & Mountain 2000). 
More recently, Davies et al. (2003) have observed 
vibration-rotation transition lines of H2 from the 
nuclei of seven ULlRGs and concluded that the 
line emission arises from dense PDRs irradiated 
by intense far-UV radiation fields. In galactic 
sources, we have used observations of H2 lines from 
high vibrational states (Luhman et al. 1995), as 
well as comparisons of extended H2 emission and 
tracers of UV excited PDRs such as far-IR contin- 
uum and [C 11] emission to show that much of the 
H2 emission arises from UV-excited gas (Luhman 
& Jaffe 1996; Luhman et al. 1997; Klumpe 1999). 
This comparison with far-IR lines and continuum 
has been more difficult in other galaxies because 
of the significant differences between the [C 11] and 
far-IR beam sizes and the typical slit sizes used 
for near-IR spectroscopy of H2. In this work, we 
will use both H2 line ratios and a study of vari- 
ous tracers of star forming clouds in the context 
of PDR models to learn more about the relative 
significance of various types of H2 excitation. 

In Section 2, we present a description of our 
spectroscopic observations of the starburst galax- 
ies M82 and NGC 253, and the normal spiral 
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galaxies NGC 6946 and IC 342. In Section 3 
we first show the observed H2 spectra and their 
extinction-corrected Une fluxes. We then show the 
spatial distribution of the H2 1-0 5(1) line inten- 
sity and compare it to the distribution of other 
PDR tracers, e.g., CO line emission and sub-mm 
continuum. In Section 4 we analyze the nature 
of the H2 excitation mechanism by comparing the 
measured H2 line ratios with thermal and non- 
thermal models. We also calculate the luminos- 
ity ratios of H2 to far-IR and [C 11] to far-IR, and 
discuss the implications of the observed ratios by 
comparison to predictions by PDR models. Sec- 
tion 5 summarizes our conclusions. 

2. OBSERVATIONS 

2.1. Fabry-Perot Spectrometer 

We have observed emission lines resulting from 
vibration-rotation transitions of H2 from the cen- 
tral 0.3 - 1.5 kpc regions of M82, NGC 253, NGC 
6946, and IC 342. Table 1 lists the basic pa- 
rameters of the four galaxies in our sample. We 
made the observations at the 1.5 m telescope of the 
Cerro Tololo Inter-American Observatory (CTIO) 
and the 5 m Hale telescope of the Palomar Ob- 
servatory^ with the University of Texas Near-IR 
Fabry-Perot Spectrometer (UTFPS). The spec- 
trometer was designed to observe extended and 
low surface brightness line emission by maximiz- 
ing the area-solid angle product (Luhman et al. 
1995). In order to subtract the telluric OH line 
emission, we chopped the secondary mirror from 
the source position to Aa ~ ±16' (at CTIO) or 
±2' (at Palomar) at 0.5 Hz. The Fabry-Perot 
etalon was scanned to cover the emission profiles, 
and each spectral step consisted of one chopping 
cycle: object-sky-sky-object (see Table 2). The 
detailed procedure for the UTFPS observations is 
described in Pak et al. (1998). The instrument, 
as used for the observations described in this pa- 
per, diff'ered in one respect from the instrument 
described in Luhman et al. (1995). For the 1996 
September run, the Fabry-Perot etalon was housed 
in an enclosure cooled thermoelectrically to — 10°C 
to reduce the background signal and to stabilize 

^Observations at Palomar Observatory were made as part 
of a continuing collaborative agreement between the Cali- 
fornia Institute of Technology, Cornell University, and the 
Jet Propulsion Laboratory 



the instrument against drifts in the plate spacing. 

The beam shape on all runs was approximately 
cylindrical (Pak et al. 1998). We indicate the 
beam size at each telescope using the equivalent 
disk size, 6ed, the diameter of a cylindrical beam 
whose solid angle is the same as the integrated 
solid angle of the actual beam pattern. Since the 
one-dimensional beam pattern was close to that of 
a box function, the FWHM of the beam pattern is 
almost the same as Oed- After the 1995 Decem- 
ber run at the Palomar Observatory, we changed 
the optical configuration inside the instrument. 
The equivalent disk size on the 1996 September 
run {6ed = H") was smaller than that on the 
1995 December run {9ed = 16"). The telescope 
pointing uncertainty was, in all cases, less than 10 
percent of the beam size. As the area-solid angle 
product of the instrument changes, the resolving 
power of the spectrometer also changes by a small 
amount, i.e., R=A/AA varies between 2300 and 
3000. Table 2 therefore also lists the effective spec- 
tral resolution for each wavelength and observing 
run. The flux calibration was done by measuring 
standard stars at the center of the beam profile. 
The flux calibration uncertainty is less than 20 
percent. Table 2 summarizes the observation logs, 
and lists the instrument parameters and the flux 
calibration stars. 

On the 1995 December run and the 1996 
September run, we used the Hale telescope, where 
the UTFPS beam has a Hnear size of 100 - 300 pc 
in our target galaxies, to map the central regions 
of these galaxies (except NGC 6946, where only 
a single point was observed). We observed H2 
1 S'(l) line emission along the major axis or the 
molecular bar known previously from CO obser- 
vations (see Table 1) stepping either 10" or 15" 
between the individual measurements. At some 
positions, we also observed 2 1 S{1) line emission. 
On the 1995 October CTIO run, we observed the 
1-0 5(1), 2-1 5(1), and 3 2 5(3) lines in NGC 253 
with a beam size of 81" (~1 kpc in this galaxy). 
This large beam covers the entire region scanned 
during the 11" beam observations in NGC 253. 

2.2. Long-slit Spectrometer 

We observed H2 emission lines from M82 with 
the IR Cryogenic Spectrometer (CRSP; Joyce et 
al. 1994) at the 2.1 m telescope of the Kitt Peak 
National Observatory (KPNO) in 1996 December. 
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The CRSP is a long-slit spectrometer with a slit 
length of 82'.'4 (135 pixels x Q/.'Ql). We used a 
300 l/mm grating and a slit width of 2.7 pixels 
(or l'/65) and the resulting spectral resolution was 
AVpwHM = 283 km s~^ at 2.2/zm, corresponding 
to R~1060. The grating was tilted to take spectra 
in the 2.11 — 2.25 /im band which includes the H2 

1- S{1) 2.12183 Atm, 1-0 S{0) 2.22329 fxm, and 

2- 1 5(1) 2.24771 /xm emission lines. 

The CRSP slit used toward M82 was rotated 
to match the major axis (P.A. = 73°) of the CO 
molecular bar (Shen & Lo 1995). In order to cover 
the extended emission from the central region, wo 
also took spectra with the slit displaced by ±2 " 
perpendicular to its length: (Aa, AS) = (— G'.'58, 
-1'.'91), (0, 0), and (+0'.'58, +1'.'91). We co-added 
9 pixels along the spatial direction (or slit length 
direction) and co-added the data taken at all three 
slit positions to make high signal-to-noise spectra 
reflecting the near-IR emission from a rectangular 
beam of 5'.'6 x 5'.'5. In the central region of M82, 
the H2 emission may extend up to 10" in the mi- 
nor axis direction (see Section 3.4.1). Even though 
the co-added rectangular beam did not include all 
of the flux from M82, the H2 line ratios from these 
data are more reliable than those derived from 
Fabry-Perot observations made at different times 
because the ratios are less dependent on the sky 
conditions (see Section 4.1). The CRSP measure- 
ments can therefore contribute to our analysis of 
the II2 excitation mechanism. 

3. RESULTS 
3.1. H2 Spectra 

Figures 1 — 7 show the H2 spectra taken with the 
UTFPS toward M82, NGC 253, NGC 6946, and 
IC 342. On each H2 spectrum, we have overlaid 
a ^^CO J = 1 — >■ spectrum taken with a similar 
beam size and smoothed in the spectral direction 
with a Gaussian smoothing function to the reso- 
lution of the UTFPS to allow us to compare the 
radial velocities and the line shapes more carefully. 
In general, those H2 line profiles with high signal- 
to-noise are similar to the CO profiles from the 
same positions. The similarity argues for a simi- 
lar distribution of line cmissivity with galactocen- 
tric radius, implying that the H2 emission and CO 
emission arise from the same molecular clouds. At 
the M82 positions (E27",N13") and (W27", S13"), 



however, the CO and H2 profiles differ. At least 
part of the change in shape and velocity profile 
may result from the beam pattern difference be- 
tween the CO and H2 observations. The spatial 
16" beam of the Nobeyama Radio Telescope is 
Gaussian while that of the UTFPS 16" beam is 
cylindrical. This shape difference could work in 
concert with the large gradient in rotational ve- 
locity at these positions to produce the profile dif- 
ferences. 

Figure 8 shows selected spectra taken by CRSP 
toward M82. The H2 1-0 5(1) lines are prominent 
and detected at all positions. The 2-1 5(1) and 

1- 5(0) lines are clearly detected at the (E10'.'5, 
N3'.'2), (W10'.'5, N3'.'2), and (W21"0, S6"4) posi- 
tions, but are confused with stellar absorption fea- 
tures at the (0, 0) position where the continuum 
level is highest. The relative strength of 1-0 5(0) 
and 2-1 5(1) in this work agrees with that of 
Forster-Schreiber et al. (2001) (see section 4.1). 
There is a possible detection of the 3-2 5(3) line 
at (W21"0, S6"4), but this line was not detected 
at other positions. 

In Figure 2, we overlaid the H2 1-0 5(l)and 

2- 1 5(l)lines from M82, as observed with the 
UTFPS. The H2 2-1 5(1) emission at (0, 0) is 
very similar to the line profile of the 1-0 5(1) tran- 
sition while the match in line shapes is less clear 
at (E14E",N6") where the signal to noise ratio of 
the 2-1 5(1) observation is small. Fluxes and up- 
per limits for the H2 2-1 5(1) and 3-2 5(3) lines 
toward NGC 253 and NGC 6946 arc listed in Ta- 
ble 3. Table 4 lists line ratios, principally for the 
2-1 5(1) to 1 5(1) ratio, along the major axis 
of M82. S/N for the ratio determination ranges 
from 2 to 6. However, temporal variations in tel- 
luric absorption may increase the uncertainty for 
the line ratios derived from the UTFPS data. 

3.2. Velocity Integrated H2 Flux 

3.2.1. Fabry-Perot Spectrometer 

We determine the H2 line flux by integrating 
across the velocity interval containing the line. We 
define the line emission bandwidth, AVj,, based on 
the 12(30 J = 1 ^ or J = 2 ^ 1 line spectra 
and list the values in Table 3. For each profile, we 
integrate over this velocity interval, subtracting a 
continuum value determined from the average of 
the channels at the two ends of the spectrum. We 
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label the flux in the portion of the spectrum where 
we are integrating over the Une as (/l.i) and that 
over the spectral region containing only contin- 
uum as {,fc,i) where the subscript i is an index for 
the individual samples along the spectrum. The 
integrated line flux (/) is calculated as: 

f = ^VchY,{fL,i-fc) , (1) 

i 

where fc is the average of fc,i and AVch is the in- 
terval between the channels in the scanning mode 
(see Table 2). The error (or noise) of the inte- 
grated line flux is derived from the standard devi- 
ation, a{fc,i), of fc,i-- 

a{f) = AVch(j{.fc,i)VNL{l + 1/Nc) , (2) 

where Ni^ and Nc are the numbers of individual 
spectral samples in the emission and continuum 
bands, respectively. Table 3 lists the observed po- 
sitions and the measured fluxes of the H2 1-0 S{1), 
2-1 5(1), and 3-2 S{3) emission lines. Note that 
the flux values, at the same wavelength and the 
same position, vary with the beam size because 
the emission regions are extended. 

3.2.2. Long-slit Spectrometer 

The FWHM {AVfwhm = 283 km s^^) of the 
instrumental profile of CRSP is wider than the 
typical FWHM of M82 emission lines. Assuming 
the line profile is gaussian, we derive the H2 line 
fiux by using the single profile fitting task in the 
IRAF^ package, SPLOT. The error in the flux is 
derived using equation (2) with 1/Nc ~ and 
AVl = 2 X AVfwhm- 

3.3. Extinction Correction 

The molecular hydrogen flux that we observe 
likely arises from the surfaces of giant molecular 
clouds exposed to far-UV radiation (Pak, Jaffe, & 
Keller 1996a, b; Davies et al. 2003). Therefore, as 
for the large-scale CO line emission from galaxies, 
the area fllling factor for the observed H2 line emis- 
sion is likely to be less than unity, so that much of 
the emergent H2 will suffer only modest extinction 

^IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 



and the foreground screen model for extinction is 
appropriate. If the H2 emission is widespread and 
arises from the surfaces of clouds, then there would 
only be a loss of half the H2 flux (from the back 
side of each cloud), and a further loss from the 
foreground extinction. The screen model is highly 
effective in modeling the extended H2 line emis- 
sion in the central kiloparsec of the Milky Way. 
Here the H2 emission would not be detectable at 
all if the emitting sources were imiformly mixed 
with molecular gas clouds, since CO observations 
indicate these clouds typically have Ay — 1500 
(Pak, Jaffe, & Keller 1996a,b). 

The far-UV radiation that gives rise to the ob- 
served H2 line emission arises from a slightly older 
population of stars than those that excite H 11 re- 
gions. These older stars are not likely to be still 
embedded in their natal clouds and can therefore 
produce the observed widespread H2 emission. In 
contrast, the Hi recombination lines likely emerge 
from Hii regions excited by very early type stars 
that are still embedded in their natal molecular 
clouds. Therefore, one might expect the extinc- 
tion of the Hi recombination lines to be signifi- 
cantly higher, and the mixed extinction scenario 
may be more likely to apply to these lines. How- 
ever, the application of a screen model even to 
the observed Br7 flux from both ultraluminous 
IR galaxies (ULIRGs) (Goldader et al. 1997b) and 
normal and nearby starburst galaxies (Kennicutt 
1998) results in a reasonably tight Br7 to far-IR 
correlation that suggests no large extinction cor- 
rections are required. 

The foreground screen model is applicable for 
all of the galaxies in our sample. The most contro- 
versial galaxy in our sample is the nearby, highly 
inclined {i — 82°) galaxy, M82 (Lynds & Sandage 
1963). Extensive modeling of the observed H re- 
combination line emission from M82 suggests Ay 
= 3 to 9 magnitudes when a uniform foreground 
screen model is applied (Lester et al. 1990; Pux- 
ley 1991; Satyapal el al. 1997; Forster-Schreiber et 
al. 2001). However, a model where the dust and 
H II regions are homogeneously distributed also flts 
the data and suggests that Ay — 25 to 52 magni- 
tudes (Puxley 1991). We argue that while a mixed 
model may apply for the Hi recombination lines, 
it is not appropriate for the H2 lines due to the 
argument outlined above. To compare the extinc- 
tion corrected results, in Table 5, we applied both 
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foreground model and mixed model to the H2 data 
of M82. 

In converting from visual extinction, Ay, in 
Table 5 to IR extinction, A\, at different wave- 
lengths, we used the IR interstellar reddening 
law (Section 7.8. in Cox (1999) and references 
therein) : 

Ax = 0.412 Av A"^-^^ (valid for 0.9 < A < 6 /xm), 

(3) 

where A is in units of jim. In section 4.1 and Ta- 
ble 4, where we compare the H2 line ratios, we 
corrected for differential extinction for lines at dif- 
ferent wavelengths. 

3.4. Intensity Distribution 

Our new H2 measurements of galaxies are par- 
ticularly useful for comparison with other obser- 
vations of line and continuum radiation from neu- 
tral atomic and molecular gas because the size of 
our UTFPS beam is more comparable to the sizes 
of beams used to measure the other PDR trac- 
ers than are the arcsec-scale beam sizes of typical 
near-IR slit spectrometers. Here we compare our 
H2 results to far-IR and submillimeter continuum, 
and [Cii] 158 /xm and ^^CO line data assembled 
from the literature. 

3.4.1. M82 

Figure 9 shows the integrated intensities of H2 
1-0 5(1), I^CO J = 1 ^ (Nakai et al. 1987), 
and i^CO J = 4 ^ 3 Hncs (Mao ct al. 2000) and 
the specific intensity of submillimeter continuum 
at 450 /zm (Hughes, Gear, & Robson 1994) along 
the major axis of M82 (P. A. = 65°; sec Tables 1 
and 3). The 450 /xm data were obtained with a 
9" beam. We smoothed these data to match the 
spatial resolution of our 16" beam observations. 
The beam sizes of the CO J = 1 — > (Nakai et al. 
1987) and the CO J = 4 ^ 3 (Mao et al. 2000) 
data arc slightly larger (16" and 18", respectively) 
than the H2 beam. 

In general, the distribution of CO, submil- 
limeter continuum, and H2 emission agree rather 
well. The differences are consistent with a com- 
mon emission region with excitation differences 
accounting for the small differences in source size. 
CO J = 1 ^ observations show emission peaks 
on either side of the nucleus (separation ~ 26"), 
indicating a rotating molecular ring (Shen & Lo 



1995). Mao et al. (2000) observed a smaller sep- 
aration (~ 15" or 270 pc) of the corresponding 
peaks in CO J = 4 — » 3. Thus, higher gas temper- 
ature and higher density clouds may exist at the 
inner edges of the rotating molecular ring. While 
the convolved scan map at 450 ixm that we show in 
Figure 9 has a plateau at the nucleus, the original 
diffraction-limited map shows weak double peaks 
separated by ~ 15" (Hughes, Gear, & Robson 
1994). This small separation supports the idea 
that the rotating molecular ring has a tempera- 
ture gradient, since the dust continuum emission 
scales at least linearly with temperature (Sohn ct 
al. 2001). The similarities in the CO, submillime- 
ter continuum, and H2 maps are consistent with a 
picture in which the H2 emission from M82 arises 
primarily from the molecular clouds in the circum- 
nuclear ring. 

3.12. NGC 253 

Figure 10 shows the integrated intensities of H2 
1-0 5(1), i^co J = 1 ^ (Canzian, Mundy, 
& ScoviUc 1988), and 12 CO J = 3 ^ 2 (Israel, 
White, & Baas 1995) as well as the specific inten- 
sity of submillimeter continuum at 850 /um (Al- 
ton, Davies, & Bianchi 1999) along the molecular 
bar in NGC 253 which has an angular extent of 
39" X 12" at a position angle of 64°. The molec- 
ular bar is tilted from the major axis (P. A. = 
51°, Pence 1980) and apparently rotates as a rigid 
body (Canzian, Mundy, & Scoville 1988). The CO 
J = 1 ^ data were obtained at the Owens Val- 
ley Radio Observatory with a 5" x 9" beam. We 
smoothed them to match our UTFPS spatial res- 
olution. The beam sizes of the CO ,7 = 3^2 line 
and the 850 jim continuum observations are 14" 
and 15", respectively. 

The H2 emission distribution in Figure 10 has 
a sharp peak at the center and extends with rel- 
atively low intensities out to ~ 0.7 kpc from the 
nucleus. The slight enhancement of II2 intensity 
at the southwestern bar (at the positions between 
—300 pc and —100 pc shown on the horizontal axis 
in Figure 10) may be related to the southwest peak 
that is also present in the overlaid CO J = 1 ^ 
distribution. The concentration of the H2 emission 
in this galaxy can be estimated by comparing the 
flux measurements with different aperture sizes. 
Engelbracht et al. (1998) derived an H2 1-0 5(1) 
flux from their 2'.'4x 12" slit of (1.21±0.07) x IQ-^'^ 
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W m ^ . From our scanned observations along the 
bar with the 11" beam (100" x 11"), the summed 
H2 1-0 5(1) flux is (2.91 ± 0.34) x lO-^^^ W m'^. 
From the 81" beam observations at the CTIO the 
flux is (3.17 ± 0.25) X 10-16 W m'^ (see Tables 
3 and 5). We conclude that the H2 emission from 
the central ^ 1 kpc region is dominated by the 39" 
long molecular bar. 

3.4.3. NGC 6946 

i^CO (J = 1 - 0) observations of NGC 6946 
show that the molecular gas is concentrated in 
the central 10" region (or 0.3 kpc at the distance 
of 5.5 Mpc, Ishizuki et al. 1990a). Using near-IR 
long-slit spectroscopy, Engelbracht et al. (1996) 
also found that the H2 emission extends up to 12" 
along the major axis. The flux in our 11" beam is 
2.4 X 10"^'' W m~^. We expect that our beam cov- 
ers most of the H2 emission in the central region 
of NGC 6946. 

3.4.4. ic 342 

An aperture synthesis map of ^^CO J = 1 — > 

in the central starburst region of IC 342 (Ishizuki 
et al. 1990b) shows a 7" (or 60 pc at the distance 
of 1.8 Mpc) diameter ring-like structure within 
the extended 40") molecular bar. Our beam 
size, 16" covers the whole molecular ring struc- 
ture. We also observed regions offset by 16" along 
the molecular bar (see Table 3). Figure 11 shows 
scan maps of the integrated intensities of the H2 
1 5(1), i^CQ J = 1 ^ (Lo et al. 1984), and 
i^CO J = 3 ^ 2 (Steppe et al. 1990) lines. Since 
the spatial resolution of the CO J = 1 ^ ob- 
servations is 4" and that of the CO J = 3 — > 2 
observations is 7", we smoothed the CO maps to 
match to our UTFPS spatial resolution of 16". 

In this galaxy, the H2 emission region is un- 
resolved and therefore significantly more compact 
than the CO emission region. This difference be- 
tween IC 342 and M82 and NGC 253 may, how- 
ever, be due to a more extreme difference in emis- 
sivities rather than to a separation of emission re- 
gions. Boker, Forster-Schreiber, & Genzel (1997) 
observed the central region of IC 342 and found 
an H2 1-0 5(1) emission region corresponding to 
the 7" CO ring. 



3.5. Total Flux 

If H2 1 5(1) line emission traces molecular 
gas illuminated by far-UV photons, its distribu- 
tion and strength should be related to those of 
other PDR probes. We have compiled far-IR con- 
tinuum, [Oi] line, and [Cii] line emission results 
from the literature. Far-IR continuum, [O i] , and 
[C 11] emission have been taken with low spatial 
resolution (0.'5 — 2'). The individual beams there- 
fore cover ~ 1 kpc regions in the centers of nearby 
galaxies (d ~ 4 Mpc;). Our big beam Fabry-Perot 
data cover comparable size scales. 

3.5.1. H2 1-0 5(1) 

The CO interferometer map of M82 (Shen & 
Lo 1995) shows that the minor-axis width of the 
rotating ring-like bar is ~ 10". Our scan map 
along the major axis with a 16" beam therefore 
most likely covers most of the H2 emission re- 
gion in the center. The total H2 1-0 5(1) flux, 
Fh2i from the 75" x 16" region within our map is 
(4.95 ±0.36) X 10-16 W m^^ (see Table 5). NGC 
253 is also an edge-on galaxy with the inclination 
of 78° (Pence 1980). The minor-axis width of the 
molecular bar, from the CO observations, is ~ 12" 
(Canzian, Mundy, & Scoville 1988), which is about 
the same as our 11" beam. For this galaxy, we 
have the 81" CTIO measurement to provide the 
total H2 flux. The active star forming regions in 
NGC 6946 and IC 342 are relatively small com- 
pared to those in the starburst galaxies, M82 and 
NGC 253. One UTFPS beam covers the H2 emis- 
sion region. 

We converted the total H2 1-0 5(1) flux values 
into extinction corrected flux values assuming that 
the foreground model is applicable (see Section 3.3 
and Table 5). 

3.5.2. Far-IR 

When making use of IRAS 60 /im and 100 /im 
observations, one can define the total far-IR con- 
tinuum flux, FpiR, as the flux measured from the 
source within a bandpass of 80 /im centered at 82.5 
lj,m (Helou et al. 1988). In addition to the IRAS 
observations, the far-IR fluxes of the galaxies in 
our II2 sample were measured using the University 
of Texas far-IR photometer (see the references in 
Smith & Harvey 1996) on the Kuiper Airborne 
Observatory (KAO) and using the Long Wave- 
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length Spectrometer (LWS; Clcgg ct al. 1996) on 
the Infrared Space Observatory (ISO; Kessler et 
al. 1996). The KAO observations spatially re- 
solved the emission regions with a beam size of 
30" X 40". The spatial resolution of the IRAS 
data, which were processed using a maximum cor- 
relation method algorithm (Aumann, Fowler, & 
Melnyk 1990), is 100" x 140" at 100 /zm, and that 
of the ISO/LWS is 80" (Gry et al. 2001). 

In Section 3.4, wc found that most of the sub- 
millimeter dust continuum emission arises from re- 
gions < 1' in diameter. Even though the mea- 
sured luminosities from the KAO instrument are 
slightly less than those from the IRAS and ISO, 
the different aperture sizes do not significantly af- 
fect the measured himinosity. Compiling the far- 
IR luminosities from the IRAS, the KAO, and the 
ISO/LWS, we list in Table 6 the average Fpm of 
the published values. 

3.5.3. [Oi] 63 jjLm and [CiiJ 158 /xm 

The [O i] 63 //m flux, Fqi, and the [C ii] 158 //m 
flux, Fcii wore measured with the U. C. Bcirkc;- 
ley cryogenic tandem Fabry-Perot spectrometer on 
the KAO (see the references in Stacey et al. 1991) 
and by the ISO/LWS. The FWHM beam sizes 
are 55" and 69" for the KAO and ISO respec- 
tively. Since the emission regions are concentrated 
in the central < 1' of the galaxies (Sections 3.4 and 
3.5.2), the flux values measured by both KAO and 
ISO can be compared with our H2 measurement of 
the galaxies. Table 6 lists the [Cii] and [Oi] line 
fluxes. 

4. DISCUSSION 

A resolution of the long-running debate about 
the excitation mechanism for the vibration-rotation 
H2 lines in galaxies and about the physical nature 
of the line formation region is critical if we are 
to reap any benefit from the diagnostic power of 
these lines. We need to examine the observational 
evidence in the context of models for the forma- 
tion and excitation of H2 to decide whether UV 
excitation followed by a radiative cascade or col- 
lisional excitation account for the population of 
the upper states of the near-IR transitions. We 
need to determine if any collisionally excited gas 
we do sec was heated by UV light incident on the 
neutral cloud or by mechanical energy in the form 



of shocks. 

Like giant molecular clouds in our own galaxy, 
the ISM of other spirals and of starburst galax- 
ies is bathed in UV radiation from main-sequence 
OB stars that have already freed themselves from 
the dense cores in which they formed. Orion 
serves as a good example for the galactic molecu- 
lar clouds. The H2 1-0 S{1) emission attributable 
to the PDRs produced by UV light has a lumi- 
nosity of 34 Lq over a 15 pc diameter region cen- 
tered on Ori. The luminosity of the same line 
attributable to the well-known and very promi- 
nent Orion/KL shock is only 2.5 Lq (Luhman et 
al. 1994). Observations of several H2 vibration- 
rotation transitions with the UTFPS, including 
detection of the 5-3 0(3) line, and comparisons 
to other PDR tracers led Pak et al. (1998) to con- 
clude that the majority of the H2 emission from 
several giant molecular clouds in the Large Mag- 
ellanic Cloud is a result of UV excitation. 

The ubiquity of UV emission incident on the 
surfaces of molecular clouds in the inner regions 
of galaxies means that these galaxies will always 
produce some UV-excited II2 emission. The real 
question then is whether this emission - or emis- 
sion from shocks produced by some combination of 
protostellar jets, supernovae, and large-scale gas 
flows - will dominate the observed flux. In this 
section, we examine whether UV excitation is suf- 
ficient to account for all of the H2 emission ob- 
served toward the galaxies in our sample. 

4.1. The H2 Line Ratios 

The ratios of H2 lines arising in states with dif- 
ferent upper state energies and the H2 ortho- to 
para-ratios give us a way to estimate the fraction 
of the emission produced by each excitation pro- 
cess within the observed beam. In PDRs where 
n(i?2) < 5 X 10^ cm~'^ (sources where a radia- 
tive cascade determines the level populations) , the 
ratios of 2-1 S{1) to 1-0 S{1) and 3-2 S'(3) to 
1-0 S{1) are ~ 0.6 and ~ 0.2, respectively (Black 
& van Dishoeck 1987). In shocked regions, e.g., 
bipolar outflows (Burton et al. 1989; Davis et al. 
1999) and supernova remnants (Richter, Graham, 
& Wright 1995), the observed ratio of 2-1 S{1) 
to 1-0 5'(1) line flux is ^ 0.1. which can be ex- 
plained by a combination of J- and C-type bow 
shocks (Smith 1991, 1994). If the level popula- 
tions of the H2 molecules were in LTE, this ratio 
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would be appropriate for gas with T^xc ~ 2000 K 
(Shull & Beckwith 1982; Black & van Dishoeck 
1987). At this temperature, the ratio of 3-2 5(3) 
to 1-0 5(1) is only ~ 0.006 (see Table 4). 

In Table 4, the observed H2 ratios of 2-1 5(1) 
to 1-0 5(1) from M82, NGC 253, and NGC 6949 
are 0.3 — 0.6, which lies between the ratios in the 
thermal and non-thermal excitation cases. If we 
assume that H2 emission arises from a mixture of 
thermal excitation at one LTE temperature and 
non-thermal excitation following the rules for a 
UV-excited cascade, we can estimate the fraction 
of the total emission in H2 lines contributed by 
UV-excited gas. In the discussion that follows, we 
use SR to represent the fraction of the H2 emission 
(summed over all H2 transitions) contributed by 
UV fluorescence. Appendix A.l gives a detailed 
explanation of how we derive from the observed 
line ratios. 

Figures 12, 13, and 14 show the values de- 
rived from our UTFPS and CRSP data from M82 
and NGC 253. We also plot the values of 3? de- 
rived from other observations of H2 emission in 
those galaxies (Forster-Schreiber et al. 2001; Har- 
rison et al. 1998; Engelbracht et al. 1998). The 
results show that, at almost all positions in both 
galaxies, the majority of the total H2 emission 
comes from regions excited radiatively by UV pho- 
tons. Figure 14 shows the values of derived from 
a variety of H2 vibration-rotation transitions for 
NGC 253. Most lines give values of 5R between 
0.6 and 1, again indicating the predominance of 
UV-excited emission. If some of the remaining 
(thermal) emission arises from dense gas in PDRs, 
the non-thermal fraction may be even higher than 
the results from the analysis used in the figures 
indicates. Sternberg & Dalgarno (1989) modeled 
the H2 emission from warm PDRs, whcirc; the H2 
vibrational level populations are thermalized by 
collisional de-excitation. In these regions, the H2 
line has excitation temperature of ^ 1000 K. If 
we use the relative line intensities derived assum- 
ing an excitation temperature, Tg^c = 1000 K for 
the collisionally excited gas, the derived K values 
are closer to unity than those derived assuming 
Te^c = 2000 K. 

Our analysis shows that II2 emission from the 
central regions of galaxies is mostly dominated by 
non-thermal excitation. UV radiation may also 
be responsible for some of the remaining thermal 



II2 emission if that emission arises in dense PDRs. 
The analysis presented in Appendix A. 2 and illus- 
trated in Figure 15 points out that the picture one 
gets from the 10 5(1) line alone is highly biased 
toward collisionally excited emission, because this 
line's upper state energy is lower than those of 
2 1 5(1) or 3 2 5(3) lines. A UV-excited radia- 
tive cascade populates many rotational states in 
higher vibrational levels. In purely radiatively ex- 
cited II2, the 1- 5(1) line accounts for only 1.6% 
of the total H2 emission, while the same line from 
2000 K gas in LTE accounts for 8.5% of the ro- 
vibrational intensity (Black & van Dishoeck 1987). 
As a result, even when only a small fraction of 
the 1-0 5(1) emission arises in UV-excited regions, 
most of the overall H2 emission is a product of UV 
excitation. In such cases, the thermal contribution 
to the 1-0 5(1) emission could arise in dense PDR 
gas where UV radiation is the ultimate source of 
the emission, or it could come from shocks associ- 
ated with supernova remnants or large-scale flows. 
A comparison of H2 1-0 5(1) luminosity with fine 
structure line emission from PDRs in galaxies will 
help distinguish between these possibilities. 

4.2. H2, [Cii], and far-IR and the Origin 
of H2 Emission 

In starburst nuclei, the far-IR luminosity arises 

mostly from PDRs illuminated by recently formed 
OB stars. Far-IR emission has been used to es- 
timate the star formation rates (SFRs) of such 
galaxies (Kennicutt 1998). The ^^3/2 p^^^ (j+ 
emission at 158 /xm from galactic nuclei also arises 
mostly from the warm photodissociated surfaces of 
far-UV exposed molecular clouds, with lesser con- 
tributions from low density ionized gas regions, 
and from the cold neutral medium (Stacey et al. 
1991). The average [C 11] to far-IR luminosity ratio 
for the four galaxies in our sample is io~2-87±o.27 
(see also Figure 16), where the 0.27 represents the 
standard deviation about the mean ratio. The 
fairly constant ratio is consistent with a PDR ori- 
gin for both (Crawford et al. 1985; Stacey et al. 
1991). 

In Figure 16, we also plot the ratio of H2 
1-0 5(1) to far-IR luminosity versus the far-IR 
luminosity from our observations. The Lh2/Lfir 
ratio is almost constant, with a mean value of 
10-5.05±o.28 Goldader et al. (1995, 1997a,b) 
found a similar proportionality between the ex- 
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tinction corrected H2 1-0 S{1) and far-IR lu- 
minosities from LIRGs (10^^ L© < Lfir < 
10^^'^ -^^©)- Their data give a value for Lh2/Lfir 
from the LIRGs of I0-5 i7±o.24^ ^j^j^j^ jg 

sistent with our value for starbursts and normal 
spiral galaxies. In calculating the mean ratio for 
nearby galaxies, we assumed that the foreground 
model was applicable when making the extinction 
correction to the M82 H2 data (see Section 3.3). 
If we applied the mixed model result, 52 magni- 
tudes of extinction, to the observed H2 1-0 S{1) 
line flux, the corresponding ratio of H2 1-0 5(1) 
line to far-IR continuum of M82 would be very 
large, Lh^/Lfir — 10~^-^ (see Tables 5 and 6). 
This ratio is two orders of magnitude larger than 
the extinction corrected line to continuum ratios 
typical for IR bright and LIRG galaxies (Goldader 
et al. 1995; Burston et al. 2001; Davies et al. 2003) 
and six times larger than ratio found for the most 
extreme H2 1-0 S{1) luminous system, NGC 6240 
(Goldader et al. 1995). 

The constant ratio of Lh2/ Lfir and Lcii / Lfir 
suggests that the far-IR, [C 11] , and H2 emission re- 
gions and their emission mechanisms are directly 
related. In examining this possibility, we assume 
that the H2 1-0 5(1) luminosity can represent the 
total H2 transition luminosity. If the II2 emission 
is mostly from PDRs and its energy level pop- 
ulations are partly thermalized, the relationship 
between the H2 1-0 5(1) and the total H2 tran- 
sition luminosity is given by equations (A8) and 
(AlO) in sections A.l and A.2. When 5ft = 0.6, 
0.8, and 1, the ratio of total H2 luminosity to 
luminosity in the I(1,0)S(1) line is 23, 34, and 
63, respectively. Since the fraction of the non- 
thermal excitation, Sft, observed in the galaxies in 
our sample is 0.8 ± 0.2, the error introduced by 
using the 1-0 5(1) line as an indicator of the total 
H2 transition intensity is a factor of 1.5 — 2 (or 
lO"'^ — 10°'^). This uncertainty is comparable to 
the standard deviation of the observed Lh2 / L fir 
ratio from normal and starburst galaxies (10°'^''), 
and LIRGs (10"-28). Using the 1 5(1) line as a 
proxy for the total H2 luminosity, we now compare 
the observed H2 emission to the intensity that in- 
dependently derived PDR models would predict 
for the same objects. 



4.3. The H2 to Far-IR Ratio Derived from 

PDR Models 

In this section, we examine the possibility that 
the observed H2 emission arises in PDRs by com- 
paring the average observed II2 1 5(1) intensity 
to the intensity predicted from PDR models with 
parameters set by observed far-IR continuum and 
[O i] and [C 11] line strengths. The emission from 
PDRs can be parameterized by the cloud density, 
riH = n-CH.) + 2n(H2), and the strength of the 
far-UV radiation field illuminating the cloud, Go 
(Tielens & HoUenbach 1985; Hollenbach & Tielens 
1999). The fine structure line ratio Fqi/Fch and 
the fine structure line to far-IR continuum ratio 
(-Fb/ + Fcii)/Ffir (see the values in Table 6) 
are relatively free from beam filling factor effects 
and can be used to derive the density and local 
UV field in the PDR. We compared the observed 
flux of [O i] , [C 11] , and far-IR (Table 6) to the PDR 
models of Kaufman et al. (1999), and list the best- 
fit values of and Go in Table 7. These values 
are different from the previous results in Stacey 
et al. (1991); Madden et al. (1993); Carral et al. 
(1994); Lord et al. (1996). One reason for the dif- 
ference is that the PDR model of Kaufman et al. 
(1999) used in this work includes additional gas 
heating due to polycyclic aromatic hydrocarbons 
(PAHs) which results in a lower predicted gas den- 
sity (Negishi et al. 2001). 

The H2 1-0 5(1) line intensity, Ih2, is also a 
function of the strength of the external far-UV ra- 
diation field. Go, impinging on the PDR, and the 
gas density, Uh, in the PDR. In order to com- 
pare this predicted H2 1-0 5(1) intensity for a 
derived nn, Go pair with the observed values, 
e.g., Lh2/Lfir in Figure 16, it is necessary to 
derive the intensity ratio of H2 1-0 5(1) to far- 
IR, Ih2/Ifir from the models. We can convert 
the derived far-UV radiation field strength, Go in 
Table 7, into far-IR intensity, Ifir, assuming that 
most of the far-UV photons are absorbed by grains 
and reradiated in the far-IR: 

Ifir = 2 x (1.3 x 10"'^) Go W m^^ sr^^ , (4) 

where the factor of 2 is for the photons with 
A > 206.6 nm. Figure 17 shows the intensity ratio, 
IH2I Ifir, as a function of Go and jih- To calcu- 
late the H2 intensity, we used a second PDR model 
that calculates abundances and temperatures in a 
plane-parallel slab for a given ur and Go (Black 
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& van Dishocck 1987; Janscn ct al. 1995) and uses 
a Monte Carlo method for the radiative transfer 
(Choi et al. 1995). We assume solar metallicity in 
the models because the H2 intensity is not sensi- 
tive to the metallicity (Pak et al. 1998). 

In Figure 17, we plot the locations of the galax- 
ies in our sample in (Go, uh) space, as derived 
from the [Oi] and [Cii] observations using the 
PDR models. The figure also shows the model 
predictions for the Ih2/Ifir ratio. Tabic 7 lists 
the values, as well as the predicted strengths of 
the H2 1-0 S{1) line. Since the derived nn val- 
ues of NGC 6946 and NGC 253 are lower than the 
model domain, we extrapolated the grids. The re- 
sulting values of Ih2/Ifir are listed in Table 7. 
The average of the H2 to far-IR ratios predicted 
by the PDR models from the fine-structure line 
observations is io-5 i3±o.36 -^^iiigi^ jg very close 
agreement with the value of io~5-05±0-28 measured 
directly (see Section 4.2). This agreement argues 
strongly for UV excitation in PDRs as the primary 
source of H2 1-0 5(1) and, by implication, as the 
dominant cause of H2 emission overall in the cen- 
tral regions of galaxies. Our results are in good 
agreement with those of Davies et al. (2003) who 
conclude that the observed H2 line emission from 
a sample of seven ULIRGs arises from dense PDRs 
exposed to strong (Go — 10^) radiation fields. 

4.4. Reasons for the Constant H2 to Far- 
IR Ratio 

In our study of the inner 400 pc of the Milky 
Way, we found that most H2 emission arises as 
a result of UV excitation and that the ratio of 
H2 1-0 S{1) luminosity to far-IR luminosity is 
- 10-4-8 (Pak, Jaffe, & Keller 1996a,b). As 
shown in sections 4.1, 4.2, and 4.3, the results 
for other galaxies are similar; UV excitation dom- 
inates and the Lh2/Lfir ratios are io-5 05±0-28 
for spiral galaxies and local starburst galaxies and 

10-5.17±0.24 foj. LIRGs. 

The relationship between far-IR intensity and 
H2 intensity as conditions in a PDR vary can be 
explained in terms of the competition between H2 
absorption and dust absorption of far-UV photons. 
Burton et al. (1990) showed analytically that the 
penetration of far-UV photons into a cloud is de- 
termined by H2 self-shielding rather than by dust 



extinction when: 

5^ > 80 cm-3. (5) 
Go 

In Figure 17, the dotted line divides the (Go, tih) 
plane into the H2 self-shielding dominant region 
and the region where shielding by dust is the ma- 
jor protector of H2. In the H2 self-shielding case, 
the UV-excitcd H2 intensity, Ih2^ is proportional 
to the far-UV field strength. Go because any addi- 
tional UV photons are absorbed by H2 and most 
of these absorptions result in a fluorescent cas- 
cade. In this case, the Ih2/Ifir ratio is always 
~ 10~^. In the case where dust shielding domi- 
nates, the H2 intensity is independent of the far- 
UV field but scales with the gas density, Uh and 
the Ih2/Ifir ratio is < 10~^. The distribution of 
our sample galaxies in the (Go, nn) plane (Fig- 
ure 17) shows that the ratio of uh/Gq, derived 
from the observations of [Cii] and [Oi] lines and 
far-IR contimmm (see section 4.3), is lower than 80 
in all cases, consistent with the observed Ih2 / Ifir 
ratio. Dust absorption of far-UV photons should 
therefore dominate over H2 self-shielding in the 
central regions of these galaxies. 

The observed constant ratio of H2 1-0 5(1) in- 
tensity to far-IR intensity is surprising in view of 
the results of the PDR modelling. In the dust- 
shielding regime, this constant ratio can only come 
about if the gas density increases as the far-UV 
field strength increases. For a uniform medium, 
with an observed \og{lH2/lFiR) = —5, the im- 
plied relationship is 

UH oc Gj-3 . (6) 

This inferred relation joins a long string of sim- 
ilar correlations between mass or density and sur- 
face density indicators and proxies for the star 
formation rate. These include: the Schmidt law 
(Schmidt 1959), as determined from Ha and CO 
J=l — » observations (Kennicutt 1998), the scal- 
ing of nu and Go derived from [O i] , [C 11] and far- 
IR intensities (Malhotra et al. 2001), the strong 
correlation of the bolometric luminosity and virial 
mass of dense galactic cores derived from CS and 
far-IR observations (Shirley et al. 2003) and the 
close correlation of bolometric luminosity with 
the amount of dense gas in galaxies (Solomon, 
Downes, & Radford 1992). One possibility, then, 
is that the relationship between density and inci- 
dent UV field inferred from the constant Ih2/Ifir 
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ratio is telling us that the surfaces of molecular 
clouds in galaxies, after stars have begun to form, 
have higher densities when the local UV field is 
higher. This conclusion is not unreasonable, given 
that less dense clumps will be ionized and photo- 
evaporated more rapidly and that the average den- 
sities of surviving clumps should be higher when 
the incident UV field is stronger. 

5. CONCLUSIONS 

Wc have analyzed near-IR H2 emission from 
the central ~ 1 kpc regions in the nearby star- 
burst galaxies, M82 and NGC 253, and the vigor- 
ously star-forming galaxies, NGC 6946 and IC 342. 
The line ratios of various H2 vibration-rotation 
lines, e.g., 2-1 ^(l), and 1-0 S'(O) to 1-0 3(1), 
show that, while both thermal and non-thermal 
excitation are responsible for the emission, non- 
thermal excitation predominates. In addition to 
the H2 line ratios, the constant ratio between H2 
1-0 S{1), [Cii], and far-IR luminosity suggests 
that the emission in the two lines and the far-IR 
continuum originate together in photodissociation 
regions. 

We also compared scan maps of H2 1-0 S{1) 
emission along the major axes (or molecular bars) 
of the galaxies with maps of CO and sub- mm emis- 
sion. The distribution of H2 line emission corre- 
lates with the sub-mm continuum emission, and 
also agrees well with the distribution of CO emis- 
sion. The H2 emission traces active star formation 
regions where the far-UV fields are intense and the 
gas is dense and warm. The correlation of H2 and 
sub-mm emission suggests that clouds have simi- 
lar structure or that the star formation and cloud 
distributions are similar in the central regions of 
these galaxies. 

We compare the luminosities of II2 1-0 S{1), 
[Cii] 158 /im, and far-IR continuum. The 
Lh2/ LpiR ratios of normal spiral and starburst 
galaxies are the same. Comparing the observed 
ratio to PDR models, wc conclude that the col- 
umn of warm H2 line emitting gas is limited by 
the extinction of far-UV photons by dust, rather 
than by II2 self-shielding in the PDRs. In this 
case, the H2 intensity is only sensitive to the gas 
density. The constant Lh2/Lfir ratio over the 
large far-IR luminosity ranges of the galaxies, then 
requires that the gas density in the PDRs increase 



as the far-UV (or far-IR) intensity increases. 
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hospitality of the staff of the Institute of Space 
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A. Derivation of Total H2 Emission 

A.l. Fraction of Non- Thermal Excitation of H2 

Assuming that H2 emission arises from a mixture of thermal and non-thermal excitation, the total H2 
intensity summing over all vibration-rotation transitions, I, and the intensity of a single transition, ly, can 
be expressed as 

I = + = ^ , and (Al) 

V 

I, = (A2) 

where the superscripts th and nth refer to the thermal and non-thermal excitation, respectively. We can 
define the fraction of non-thermal H2 intensity in the observed H2 intensity as 

Tnth 

SR = — , and (A3) 

Tnth 

SR. = (A4) 

The fractional intensity of thermal or non-thermal emission in a single transition relative to the total thermal 
or non-thermal emission is: 

jth 

Sj; = , and (A5) 

Tnth 

= fru- (A6) 

Table 8 lists the adopted values of s„ and for H2 1-0 S{1), 2-1 S{1), 3-2 5(3), and 1-0 5(0) from model 
calculations (Black & van Dishoeck 1987). We assumed that the excitation temperature is 2000 K for 
thermalized gas. From the above equations we can derive the H2 line ratio as 

Iy_ ^ Sy - (.s„ -p,,)Sft 

where the subscript '1' denotes the 1-0 5(1) transition line. From equation (A7), we can rewrite, 5ft, as a 
function of the observed H2 ratio of another transition to the 1-0 5(1) line: 

[Sv -Pv) - (si -pi)[Iv/h) 

A.2. Fraction of H2 10 5(1) Line Intensity 

The H2 1-0 5(1) at 2. 1218 /um line is the most commonly observed near-IR H2 line because the wavelength 
resides in a clear part of the atmospheric window and the line is usually brighter than other H2 lines. When 
the level populations of excited H2 molecules arc partly thermalized, the 1-0 5(1) emission is enhanced 
relative to most other H2 lines because of its lower excitation requirements. We need to check whether the 
H2 1-0 5(1) emission can represent the total H2 intensity. 

We can generalize the fraction of the non- thermal excitation in a single transition, Sft„, in terms of 3? as: 



Pv V 



(A9) 
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Figure 15 shows SR^, versus K for various transitions. For example, when 5R 2± 0.6, about 20 percent of the H2 
1-0 S{1) emission arises from non-thermal excitation. On the other hand, we can extrapolate to the total 
H2 vibration-rotation transition intensity, /, from a single H2 transition line emission, ly , as 
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Fig. 1.— M82 spectra of the H2 1-0 S{1) line ob- 
served with the UTFPS of a 16" beam (soUd lines). 
The dotted lines show ^^CO J = 1 ^ spectra 
(Nakai et al. 1987) taken with the 16" beam of the 
Nobeyama 45 m telescope and smoothed to the 
spectral resolution of the H2 spectra. The R.A. 
and Dec. offsets in the upper left corner of each 
spectrum are in arc seconds from the center of the 
galaxy (see Tables 1 and 3). The Y axis gives 
the intensity scale for the H2 spectra while the in- 
tensity of each CO spectrum has been scaled to 
match the corresponding H2 spectrum. The two 
thick vertical ticks on the X axis mark the emis- 
sion bandwidth, AVl , within which we integrated 
the emission line fluxes (see Table 3). 



Fig. 2.— M82 spectra of the H2 1-0 S{1) (left 
panel) and 2 1 S{1) (right panel) lines observed 
with the UTFPS of a 11 " beam. The dotted lines 
show CO J = 1 — > spectra taken with a 16 " 
beam (see Fig. 1). 



Fig. 3.— NGC 253 spectra of the H2 1-0 S{1) line 
observed with the UTFPS of a 11" beam (solid 
hues). We have overlaid ^^CO J = 2 ^ 1 ob- 
servations obtained by Mauersberger et al. (1996) 
(beam size 12", dotted lines), smoothed in the 
spectral direction to the resolution of the H2 spec- 
trum). See Fig. 1 for further explanation. 



Fig. 4.— NGC 253 spectra of the H2 1-0 ^(l) (left 
panel) and 2-1 S{1) (right panel) lines observed 
with the UTFPS of a 11 " beam. See Figures 1, 
2, and 3 for further explanation. 



Fig. 5.— NGC 253 spectra of H2 1-0 S{1) (left 
panel), 2-1 ^(l) (middle panel), and 3-2 5(3) 
(right panel) lines observed toward the center of 
the galaxy with the UTFPS of an 81" beam. We 
have overlaid ^^CO J = 1 ^ observations ob- 
tained by Paglione et al. (2001) with the 45" beam 
of the FCRAO 14 m telescope, smoothed to the 
spectral resolution of the H2 observations. 



Fig. 6.— NGC 6946 spectra of H2 1-0 S{1) (left 
panel) and 2-1 S{1) (right panel) lines observed 
toward the center of the galaxy with the UTFPS 
of an 11" beam. The dotted line shows ^^CO J = 
1 — > observations with the 16" beam (Sofue, Doi, 
& Ishizuki 1988). 



Fig. 7.— IC 342 spectra of H2 1-0 S{1) line (solid 
line) obaserved with the UTFPS of a 16" beam. 
We smoothed the ^^CO J = 2 ^ 1 observations 
(dotted line; Xie, Young, & Schloerb 1994), taken 
with a 23" beam, to the spectral resolution of the 
H2 results. 



Fig. 8. — M82 spectra observed with the long slit 
spectrometer, CRSP. Each spectrum shows the av- 
erage over a 5"6 x 5"5 box. The positions of the H2 
1-0 5(1), 3-2 S{3), 1-0 S{0), and 2-1 S{1) lines 
are marked by arrows. Offsets for the spectra are 
in arc seconds from the central position given in 
Table 1. 



Fig. 9. — Intensity distributions of H2 1-0 S{1) 
(solid line; this work), 450 continuum (dot- 
ted line; Hughes et al. 1994), ^^CO J = 1 ^ 
(dashed line; Nakai et al. 1987), and ^^CO J = 
4^3 (dash-dotted line; Mao et al. 2000) from 
M82 along the major axis (P. A. = 65°). Positive 
numbers along the horizontal axis increase to the 
northeast on the sky. We convolved the 450 /itm 
data to match the size of the H2 and CO beams. 
The ordinate is in units of integrated H2 intensity. 
The intensity levels of the other tracers are nor- 
malized to match the H2 peak. The error bars for 
the H2 data give the intensity uncertainty at each 
point. 



Fig. 10. — Intensity distribution along the molec- 
ular bar of NGC 253 of H2 1-0 S(l) (sohd hne; 
this work), compared to the 850 fiia continuum 
(smoothed to 16 " resolution, dotted line; Alton et 
al. 1999), i^CO J = 1 ^ (dashed line; Canzian 
et al. 1988), and ^^CO J = 3 ^ 2 (dash-dotted 
line; Israel et al. 1995). 



19 



Fig. 11. — Intensity distributions of H2 1-0 S{1) 
(solid line; this work), ^^CO J = 1 — » (dashed 
line; Lo et al. 1984), and ^^CO J = 3 ^ 2 (dash- 
dotted line; Steppe et al. 1990) along the IC 342 
molecular bar (P. A. = 14°). 



Fig. 12. — Distribution of the derived ratio of 
non-thermal to total ro-vibrational H2 omission 
(see eq. [A8]), 5R, based on the observed ratios of 
H2 2 1 S{1) to 1 5(1) (top plot) and 1 5(0) 
to 1 5(1) (bottom plot) in MS 2. The large 
filled circles denote values derived from UTFPS 
11" beam data. The small filled circles and the 
small filled squares denote derived from CRSP 
5"5 X 5'.'6 aperture data. The open diamonds 
and the open squares denote 3? derived from line 
ratios obtained with a 2'.'5 x 2'.'5 (or 44 pc x 
44 pc) aperture by Forster-Schrcibcr et al. (2001). 
The abscissa shows the projected positions from 
the center along the major axis (P.A. of 65° for 
the UTFPS data and the Forster-Schreiber et al. 
(2001) data and 73° for the CRSP data). 



Fig. 13. — Distribution of 3? derived from the 
observed ratios of H2 2 1 5(1) to 1-0 5(1) (top 
plot) and 1-0 5(0) to 1-0 5(1) (bottom plot) in 
NGC 253. The abscissa shows the projected posi- 
tions from the center along the major axis (P.A. = 
51°). The filled circle denotes derived using the 
UTFPS 11" beam data. The open diamonds and 
the open squares denote show values derived from 
spectra taken with a 3" x 3" (or 36 pc x 36 pc) 
aperture by Harrison et al. (1998). 



Fig. 14. — 5R values derived from the observed ra- 
tios of various H2 lines to 1-0 5(1) line in NGC 
253. The filled circles denote the UTFPS 81" (or 
980 pc) beam data from this work. The open 
squares denote the grating spectrometer data with 
a 2'.'4 X 12" (or 29 pc x 140 pc) aperture by En- 
gelbracht et al. (1998). All data in this plot were 
observed toward the center of the galaxy. The ab- 
scissa shows the wavelength of the H2 line that 
was compared to the 1-0 5(1) line. 



Fig. 15. — Plot of versus using equa- 
tion (A9). 3?^ is the fraction of non-thermal 
H2 emission intensity in the particular line (e.g., 
1-0 5(1), 1 5(0), 2-1 5(1), or 3-2 5(3)) indi- 
cated against the plotted curve. 



Fig. 16.— Ratios of [Cii] and H2 1-0 5(1) line 
luminosity to far-IR continuum luminosity toward 
the galaxies in our sample. The [C 11] and H2 lu- 
minosity ratios are denoted by the open diamonds 
and the filled circles, respectively. See Table 6 for 
references and a detailed explanation. 



Fig. 17.— Ratios of H2 1-0 5(1) fine intensity 
to far-IR continuum intensity as a function of the 
gas density, nn, and the far-UV strength. Go. The 
contour labels are in log scale. The dotted line di- 
vides the (Go, plane into the H2 self-shielding 
dominant region and the dust absorption domi- 
nant region (see equation [5]). The sample galax- 
ies, whose Go and Uh are derived from the [Cii] 
and [O i] observations, are overlayed. 
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Table 1 

Adapted Parameters for Galaxies 



Galaxy 


a(J2000) 


5(J2000) 


ref 


P.A. 


ref 


Distance 


ref 




m s 'j 


(0 , ») 




n 




(Mpc) 




M82 


09 55 52.6 


+69 40 47 


1 




2 


3.6 


3 


NGC 253 


00 47 33.1 


-25 17 16 


4 




5 


2.5 


6 


NOG 6946 


20 34 52.2 


+60 09 14 


7 


62^ 


8 


5.5 


9 


IG 342 


03 46 48.1 


+68 05 46 


10 


14b 


11 


1.8 


12 



Position angle of the major axis. 
''Position angle of the molecular bar. 



References. — (1) Rickc et al. 1980; (2) Nilson 1973; (3) Frccdman 1994; (4) 
Becklin, Fomalont, & Ncugcbaucr 1973; (5) Caiizian, Mundy, & Scovillc 1988; (6) 
Maucrsbcrgcr ct al. 1996 (references therein); (7) Turner & Ho 1983; (8) Rogstad 
& Shostak 1972; (9) Ishizuki ct al. 1990a (references therein); (10) Becklin et al. 
1980; (11) Ishizuki et al. 1990b; (12) Boker et al. 1999 (references therein) 



Table 2 

Observation Logs and Instrument Parameters for Fabry-Perot Spectrometer 



Galaxy 


H2 Line^ 
(/iin) 


Observation 
Date 


Observatory^ 


(") 


AVfwhm 


AVbw" 
(km s"^) 




Flux Calibration 
Star 


M82 


1-0 S(l) 2.12183 


1995 Dec 


Palomar 


16 


100 


156 


43.8 


HR 8505 






1996 Sop 


Palomar 


11 


105 


164 


35.0 


HR 7924 




2-1 S(l) 2.24771 


1996 Sop 


Palomar 


11 


ill 


173 


35.0 


HR 8781 


NGC 253 


1-0 S(l) 2.12183 


1995 Oct 


ctio 


81 


125 


192 


43.8 


HR 6084, 7525, 8728 






1996 Sop 


Palomar 


11 


105 


164 


35.0 


HR 7924 




2-1 S(l) 2.24771 


1995 Oct 


ctio 


81 


132 


203 


35.0 


HR 6084 






1996 Sep 


Palomar 


11 


111 


173 


35.0 


HR 8781 




3-2 S(3) 2.20139 


1995 Oct 


CTIO 


81 


129 


198 


43.8 


HR 6084, 8728 


NGC 6946 


1-0 S(l) 2.12183 


1996 Sop 


Palomar 


11 


105 


164 


35.0 


HR 7924 




2-1 S(l) 2.24771 


1996 Sop 


Palomar 


11 


111 


173 


35.0 


HR 8781 


IC 342 


1-0 S(l) 2.12183 


1995 Dec 


Palomar 


16 


100 


156 


25.0 


HR 8505 



^Rcst wavelengths in vaeua. 

^We used the CTIO 1.5 m telescope and the Palomar 5 m telescope. 
•^Diameter of the equivalent disk. See the text for details. 

■^Full-width at half-maximum of the instrument spectral profile for an extended source. 
■^Band width of the instrument profile. AVbw — f I{y)(iV/ Ip^ak 
*^Step distance between the wavelength channels in the scanning mode. 
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Table 3 
Observed Flux of H2 Lines 



Galaxy e 


D 


Position 




1-0 S(l) 


2-1 S(l) 3-2 S(3) 












2.12183 ixm 


2.24771 ixm 2.20139 jiun 


(") 


(pc) 


(") 




(km s"^) 




(10-" W m-^) 


M82 11.1 


194 


14 


6 


+80, +430 


4.00 ± 1.04 


2.61 + 0.73 












+80, +430 


11.79 ± 2.05 


9.28 + 0.87 


15.7 


274 


41 


19 


+ 130, +480 


-1.10 + 0.81 








27 


13 


+ 130, +480 


3.74 ± 0.39 








14 


6 


+80, +430 


13.12 ± 1.71 














+80, +430 


15.99 ± 2.51 








-14 


-6 


0, +350 


12.09 ± 1.64 








-27 


-13 


0, +350 


4.51 ± 1.03 








-41 


-19 


0, +350 


-1.10 ± 1.19 




NGC 253 11.1 


135 


45 


22 


+40, +360 


0.53 ± 1.52 








36 


18 


+40, +360 


2.81 ± 0.79 








27 


13 


+40, +360 


0.71 ± 1.09 








18 


9 


+40, +360 


1.62 ± 1.09 








9 


4 


+40, +360 


2.05 ± 0.86 














+70, +390 


9.71 ± 0.91 


3.53 ± 1.54 






-9 


-4 


+ 100, +420 


6.87 ± 0.79 








-18 


-9 


+ 100, +420 


2.69 ± 1.04 








-27 


-13 


+ 100, +420 


1.22 ± 0.89 








-36 


-18 


+ 100, +420 


0.03 ± 1.22 








-45 


-22 


+ 100, +420 


0.86 + 0.73 




81 


982 








+70, +390 


31.7 ±2.5 


20.9 + 9.2 5.8 + 2.9 


NGC 6946 11.1 


296 








-60, +250 


2.42 ± 0.53 


0.75 ± 0.40 


IC 342 15.7 


137 


5 


15 


-70, +120 


-0.11 ± 0.37 






-15 



-70, +120 
-70, +120 



1.71 ± 0.39 
0.13 ± 0.46 



^Offset from the centers listed in Table 2. 

'^Emission line range (minimum, maximum) for the integrated line flux. See equation (1). 
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Table 4 
Observed H2 Line Ratios'' 





Obd 
V ) 


("■) 


Ai5 


1-0 S(0)/l-0 S(l) 


2-1 S(l)/l-0 S(l) 


3-2 S(3)/l-0 S(l) 


Non-thermal Model'' 








0.46 




0.56 




0.18 


Thermal ModeF 








0.21 




0.082 




0.0057 


M82 


5.5 X S-e** 


31.6 


9.6 


0.21 ± 0.12 





40 ±0 


13 








26.3 


8.0 


0.23 ±0.07 





20 ±0 


04 








21.0 


6.4 


0.24 ± 0.06 





26 ± 


06 








15.8 


4.8 


0.35 ± 0.05 





29 ± 


04 








10.5 


3.2 


0.33 ± 0.05 





32 ± 


05 








5.3 


1.6 


0.44 ± 0.13 





33 ± 


06 








0.0 


0.0 


0.38 ± 0.21 





47 ± 


10 








-5.3 


-1.6 


0.39 ±0.14 





23 ±0 


06 








-10.5 


-3.2 


0.38 ±0.11 





33 ±0 


05 








-15.8 


-4.8 


0.47 ± 0.06 





30 ± 


03 








-21.0 


-6.4 


0.45 ± 0.06 





27 ± 


04 








-26.3 


-8.0 


0.65 ± 0.20 





29 ± 


11 








-31.6 


-9.6 


0.35 ±0.15 





21 ± 


09 






11.1 


14 


6 







62 ± 


24 



















75 ±0 


15 




NGC 253 


11.1 













34 ±0 


15 






81 













61 ±0 


27 


0.18 ±0.09 


NGC 6946 


11.1 













30 ±0 


17 





^We correct the extinction difference for each line, assuming the foreground extinction model (see Table 5). 
^Fluorescent H2 emission for the PDR model 14 in Black & van Dishoeck (1987). 

"^Thermally excited H2 emission in LTE with Te^^c ^ 2000 K (Shull & Beckwith 1982; Black & van Dishoeck 1987). 
■^Size of the rectangular box of b'.'b (major axis) X 5'.'6 (minor axis). 



Table 5 

Total H2 1-0 S{1) Flux and Extinction Correction 



Galaxy 


Obs. Flux^ 




rcf 


Corr. Flux 




(10-" W m-2) 


(mag) 




(10-" W 


M82 


49.5 


5.5<= 


1,2 


86.9 






52^ 


3 


9820 


NGC 253 


31.7 


7.8'= 


4 


69.9 


NGC 6946 


2.42 


4.3^= 


5 


3.74 


IC 342 


1.71 




6 


3.68 



''Total H2 10 5(1) flux from the observed regions in the galaxy. 

''V-band extinction in magnitudes. 

"^Based on the foreground extinction model. 

Based on the mixed extinction model. 

"Average value from the four regions in Boker, Forster-Schreiber, & 
Genzel (1997). 

References. — (1) Lester ct al. 1990; (2) Puxlcy 1991; (3) Forster- 
Schreiber et al. 2001; (4) Engelbracht et al. 1998; (5) Engelbracht et al. 
1996; (6) Boker, Forster-Schreiber, & Genzel 1997 
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Table 6 

Observed Flux of Far-IR, [Oi], and [Cii] 



Galaxy 


Fh2 


ref 




ref 


Fell 


ref 


Poi 


ref 




(10-" 


W m- = ) 


(10-1= 


W m-=) 


(10-1^ 


W m- = ) 


(IQ-i^ 


W m-=) 


M82 


86.9 


1 


82.1 


2,3,4,5,6,7 


128 


4 


169 


4 


NGC 253 


69.9 


1 


64.7 


4,5,6,7 


51.9 


4 


37.6 


4 


NGC 6946 


3.74 


1 


3.68 


3,4,5,6,7 


10.3 


4 


5.90 


4 


IC 342 


3.68 


1 


12.2 


2,3,5,7 


38.7 


8 


32.6 


8 



^Average over the published values. 



References. — (1) This work; (2) Crawford et al. 1985; (3) Staeey et al. 1991; (4) Negishi 
et al. 2001; (5) Young et al. 1989; (6) Smith & Harvey 1996; (7) Mouri & Taniguchi 1992; (8) 
Eckart et al. 1990 



Table 7 
Results from FDR Models 



Galaxy 


log uh ° 

(em-^*) 


log Go 


log Ih2/Ifir ^ 


M82 


2.7 


2.9 


-5.2 


NGC 253 


1.8 


2.6 


-5.6 


NGC 6946 


2.5 


2.3 


-4.9 


IC 342 


2.8 


2.5 


-4.8 



'''Deduced physical parameters by comparing to the 
model results in Kaufman et al. (1999). 

^Ratio of H2 1-0 5^(1) to far-IR intensity. Ih2 is de- 
rived from the PDR model in section 4.3 and Ifir is 
from equation (4). 



Table 8 

Portion of a H2 Line in the Total H2 Emission 



Model 


Symbol 


1-0 S(l) 


1-0 S{0) 


2-1 S(l) 


3-2 5(3) 






2.12183 //m 


2.22329 ixm. 


2.24771 ixm. 


2.20139 /im 


Non-thermal ^ 




0.016 


0.0073 


0.0089 


0.0028 


Thermal 




0.085 


0.018 


0.0070 


0.00048 



*Prom the PDR model 14 in Black & van Dishoeck (1987). 

''Prom the LTE model with Texc = 2000 K (Black & van Dishoeck 1987). 

fnth I jnth 

Pv — J-v / 

d „ T^h I jth 

Sv — /Jj. 
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